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Neutron  Diffraction  Investigations  of  Magnetism  in  BiFe03 
Epitaxial  Films 


William  Ratcliff II*  Daisuke  Kan*  Wangchun  Chen,  Shannon  Watson,  Songxue  Chi, 
Ross  Erwin,  Garry  J.  McIntyre,  Sylvia  C.  Capelli,  and  Ichiro  Takeuchi 


The  recovery  of  a  modulated  magnetic  structure  in  epitaxial  BiFeOs  thin  films 
as  revealed  by  neutron  diffraction  is  reported.  The  magnetic  structure  in  thin 
films  is  found  to  strongly  depend  on  substrate  orientation.  The  substrate 
orientation  causes  different  strain-relaxation  processes  resulting  in  dif¬ 
ferent  thin-film  crystal  structures.  The  (110)  oriented  film  with  a  monoclinic 
structural  phase  has  a  single-domain  modulated  magnetic  structure  where 
the  magnetic  moment  lies  in  the  HHL  plane.  For  the  (111)  oriented  film  that 
has  a  rhombohedral  structure,  a  modulated  structure  superimposed  on  the 
G-type  antiferromagnetic  order  is  found.  These  results  indicate  that  slight 
structural  modifications  in  the  BiFe03  thin  film  cause  drastic  changes  in  the 
magnetic  structure. 


1.  Introduction 

Multiferroics  are  materials  that  exhibit  both  magnetic  and  fer¬ 
roelectric  order  parameters.  Recent  years  have  seen  a  renais¬ 
sance  in  the  study  of  these  materials,  resulting  in  the  discovery 
of  many  new  materials  in  which  the  magnetic  and  ferroelec¬ 
tric  order  parameters  are  coupled.  However,  for  most  of  these 
materials,  the  ordering  temperatures  are  far  too  low  for  device 
applications.  Perovskite-structured  BiFe03  (BFO)  with  a  rhom¬ 
bohedral  distortion  along  the  [111]  direction  (space  group  R3c) 
is  a  rare  exception  to  this  trend,  with  a  ferroelectric  transition 
temperature  Tc  of  820  “Cl1!  and  an  antiferromagnetic  transition 
temperature  TN  of  370  °C.I2I  It  has  been  demonstrated  that  an 
applied  electric  field  can  effect  the  magnetic  domain  population 
in  BFO  thin  films!3!. 
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The  study  of  the  magnetic  structure  of 
bulk  BFO  started  with  powders  and  evolved 
to  single  crystals  when  large  ones  became 
available.  These  single-crystal  studies  defin¬ 
itively  showed  that  the  magnetic  structure 
in  the  bulk  was  a  long-wavelength  spiral!4! 
in  the  plane  which  contains  one  of  the 
ferroelectric  polarization  and  the  propa¬ 
gation  vector.  The  rhombohedral  crystal¬ 
lographic  symmetry  results  in  a  degen¬ 
eracy  in  the  free  energy,  giving  rise  to  the 
possibility  of  magnetic  domains.  Each 
possible  magnetic  domain  is  character¬ 
ized  by  one  of  three  equivalent  propaga¬ 
tion  vectors  (related  by  rotations  about 
the  3  fold  axis):  fet  =  [309] , ^2  =  [039], 
k3  =  [330] ,  with  3  =  0.0045 ,  which  corresponds  to  a  periodicity 
of  620  A.  I4!  Remarkably,  as-grown  bulk  crystals  form  with  only  a 
single  ferroelectric  and  a  single  magnetic  domain.  Furthermore, 
associated  with  each  propagation  vector,  there  is  the  possibility 
that  the  spiral  could  be  either  left  or  right  handed.  However,  the 
spirals  in  the  bulk  were  found  to  possess  a  single  chirality. I4! 

However,  there  have  not  been  similarly  detailed  studies  of  the 
magnetic  structure  in  thin  film  samples.  Such  work  is  critically 
important,  as  the  multiferroic  properties  of  BFO  are  likely  to  be 
highly  dependent  on  the  changes  in  strain  that  accompany  the 
transition  from  bulk  crystal  to  thin  film  geometry.  One  recent 
study!5!  hag  shown  that  the  spiral  as  observed  in  the  bulk  is  not 
present  in  the  thin  film.  This  means  that  the  long-wavelength 
spiral  observed  in  the  bulk  collapses  to  a  G-type  antiferromag¬ 
netic  order  in  which  the  Fe  magnetic  moments  are  coupled  fer- 
romagnetically  in  the  (111)  planes  and  antiferromagnetically 
between  adjacent  planes,  leading  to  the  doubling  of  the  mag¬ 
netic  unit  cell  in  all  crystallographic  directions  from  the  nuclear 
unit  cell.!22!  This  is  manifested  as  a  single  peak  at  the  (0.5  0.5  0.5) 
position  in  reciprocal  space.  However  the  detailed  magnetic 
structure  in  the  thin  film  form  of  BFO  still  remains  unclear. 

In  this  paper,  we  investigate  how  the  thin  film  geometry 
affects  the  magnetic  structure  of  BFO.  To  this  end,  we  pursued 
a  neutron  diffraction  study  of  the  magnetic  structure  of  BFO 
epitaxial  films  grown  on  SrTi03  (STO)  substrates  with  various 
orientations.  Our  unpolarized  and  polarized  neutron  diffraction 
studies  reveal  that  the  film  has  a  different  magnetic  structure 
depending  on  the  substrate  orientation.  For  the  (110)  oriented 
film  with  the  monoclinic  Mb  structural  phase,  the  magnetic 
structure  is  found  to  be  a  single-domain  modulated  magnetic 
structure  where  the  magnetic  moment  lies  in  the  HHL  plane. 
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The  (111)  oriented  film,  in  which  the  rhombohedral  structure 
is  maintained,  has  a  modulated  structure  superimposed  on  the 
G-type  order.  These  results  indicate  that  the  magnetic  structure 
of  BFO  thin  films  is  susceptible  to  subtle  changes  in  its  crystal 
structure  resulting  from  the  different  epitaxial  strain  relaxation 
process  dependent  upon  the  substrate  orientation. 

2.  Experimental  Results  and  Discussions 

Due  to  the  rhombohedral  distortion  along  the  [111]  direc¬ 
tion,  epitaxially-grown  BiFe03  thin  films  have  crystallographic 
domain  structures  depending  on  substrate  orientations.  The 
(001)  oriented  film,  for  instance,  has  four  equivalent  [111]  crys¬ 
tallographic  directions,  creating  four  crystallographic  domain 
structures.  However,  this  multiple  crystallographic  domain 
structure  makes  determination  of  the  magnetic  structure  com¬ 
plex,  since  the  multiple  magnetic  domain  structures  could 
form  even  in  a  single  crystallographic  domain  due  to  a  six-fold 
degeneracy  for  the  orientation  of  the  magnetic  moment  within 
a  given  (111)  plane.!3,6!  To  minimize  this  complexity,  we  have 
chosen  thin  films  grown  on  (110)  and  (111)  STO  substrates  in 
which  a  (quasi-)  single  crystallographic  domain  structure  can 
be  achieved  as  confirmed  in  Figure  1.  The  BFO  films  studied 
here  have  been  fabricated  by  pulsed  laser  deposition  with  a 
thickness  of  1  pm  in  order  to  minimize  the  strain  effect  which 
could  modify  the  magnetic  structure  in  BFO.!7!  Our  recent 
X-ray  diffraction  study!8!  confirms  that  the  fabricated  BFO  films 
have  different  structures  depending  on  the  substrate  orienta¬ 
tions  even  if  the  film  thickness  exceeds  the  critical  thickness 
of  -100  nm,l910l  at  which  BFO  thin  films  begin  to  relax  from 
the  epitaxial  strain  induced  by  the  substrate.  For  the  (llO)-ori- 
ented  film,  the  crystal  structure  is  found  to  be  a  monoclinic 
structure,  the  so-called  Mb  phase.!8,11!  The  lattice  parameters  are 
determined  to  be  am  =  5.645  A,  bm  =  5.589  A,  cm  =  3.950  A  and 
P  =  90.85°  (The  subscript  m  denotes  monoclinic).  In  the  recip¬ 
rocal-space-maps  (RSMs)  shown  in  Figure  la  and  b,  we  see 
two  reflections  corresponding  to  two  crystallographic  domains 
in  the  film.  However,  it  should  be  noted  that  one  reflection 
has  an  intensity  more  than  twenty  times  larger  than  the  other. 
(The  intensity  ratio  slightly  depends  on  the  reflection  index.) 
Therefore  the  contribution  from  the  minor  domain  is  neg¬ 
ligible  and  the  film  can  be  regarded  as  having  a  quasi  single- 
crystallographic  domain.  In  the  case  of  the  (lll)-oriented  film, 
a  single  crystallographic  domain  structure  with  rhombohedral 
symmetry  is  found  to  be  maintained  and  the  lattice  parameters 
are  determined  to  be  a  =  3.968  A  and  a  =  89.4.  These  lattice 
parameters  are  in  good  agreement  with  the  lattice  parameters 
in  bulk.!12!  Figure  lc  shows  a  RSM  around  the  STO  (132)  Bragg 
reflection  taken  for  a  (lll)-oriented  thin  film.  The  observed 
single  reflection  from  the  BFO  layer  indicates  that  the  no  twin 
crystallographic  domain  is  formed  along  the  [-110]  direction. 
We  note  that  around  the  (114)  Bragg  reflection  (along  the  [11-2] 
direction)  there  is  also  only  a  single  reflection  observed  (not 
shown).  These  observations  reveal  that  the  (111)  film  is  a  single 
crystallographic  domain.  Ferroelectric  properties  for  the  fabri¬ 
cated  BFO  films  are  confirmed  by  room-temperature  square¬ 
shaped  ferroelectric  hysteresis  loops  as  seen  in  Figure  Id. 
The  observed  saturated  ferroelectric  polarizations  are  -85  and 


Figure  1.  Structural  and  ferroelectric  properties  for  (110)  and  (111)  ori¬ 
ented  BiFe03  thin  films,  a-c)  show  semi-logarithmic  contour  plots  of  the 
X-ray  reciprocal  space  mappings  taken  for  lpm-thick  BiFe03  (BFO)  films 
grown  on  SrTi03  (STO)  substrates  with  (110)  and  (111)  orientations.  The 
mappings  in  (a)  and  (b)  were  taken  for  the  (110)  film  around  the  STO 
(222)  and  (220)  Bragg  reflections,  respectively.  Note  that  the  splitting 
in  reflection  intensity  was  observed  only  along  the  [001]  direction.  The 
two  reflections  observed  from  the  BFO  layer  indicate  that  two  crystallo¬ 
graphic  domains  exist  in  the  film.  Based  on  the  reflection  intensity  ratio, 
the  population  of  one  of  the  domains  is  negligibly  small.  In  (c),  the  scans 
were  performed  with  the  (111)  oriented  film  around  the  STO  (132)  Bragg 
reflection,  exhibiting  a  single  reflection  from  BFO  layer.  This  confirms 
that  the  (1 1 1  [-oriented  film  has  the  single  crystallographic  domain  struc¬ 
ture.  In  (a  -c),  STO,  SRO,  and  BFO  denote  SrTi03,  SrRu03,  and  BiFe03, 
respectively.  The  mappings  are  plotted  in  units  of  the  reciprocal  lattice 
unit  of  STO.  d)  Room-temperature  polarization  hysteresis  loops  for  the 
(110)  and  (111)  oriented  films.  The  loops  were  measured  at  25  kHz. 


-110  pC  cm-2  for  the  (110)  and  (lll)-oriented  film,  respectively. 
These  values  are  in  good  agreement  with  the  value  predicted  by 
the  projection  rule.!13! 

We  used  neutron  scattering  to  probe  the  magnetic  struc¬ 
tures  of  the  films.  The  equation  for  the  differential  cross  sec¬ 
tion  describing  the  coherent  elastic  scattering  of  neutrons  from 
magnetically  ordered  crystals  in  the  ground  state  is  given  by: 

[14-16] 


where  NM  is  the  number  of  magnetic  unit  cells  in  the  crystal 
and  VM  denotes  the  volume  of  the  magnetic  unit  cell.  The  mag¬ 
netic  structure  factor  FM  is  defined  by: 
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where  the  neutron-electron  coupling  constant  in  parentheses 
is  -0.27  x  10-14  m,  Gand  jyf  are  unit  vectors  in  the  direction  of 
the  reciprocal  lattice  vector  G  and  the  average  moment  direc¬ 
tion.  Here  (p})  is  the  thermal  average  of  the  ordered  magnetic 
moment  of  the/*1  atom  in  the  unit  cell,  ri  is  the  position  of  the 
j™  atom  in  the  magnetic  unit  cell,  fj(Q)  is  the  scalar  magnetic 
form  factor  of  the/*1  atom  in  the  cell,  and  the  sum  j  ranges  over 
all  atoms  in  the  unit  cell.  The  (scalar)  magnetic  form  factor  is 
the  Fourier  transform  of  the  magnetization  density  associated 
with  each  atom.  An  important  consequence  of  the  factor  in  tri¬ 
angular  brackets  in  (2)  is  that  the  observed  cross  section  is  only 
sensitive  to  the  component  of  the  moment  in  the  plane  perpen¬ 
dicular  to  the  measured  reflection. 

We  shall  start  our  discussion  with  neutron  diffraction  meas¬ 
urements  of  the  1  gm-thick  film  grown  on  a  (110)  oriented  STO 
substrate.  Due  to  the  close  lattice  constants  of  STO  (a  =  3.905A) 
and  BFO  (a  =  3.%A,  a  =  89.4°  in  buk)  and  the  relatively  weak 
signal  from  the  film  relative  to  the  substrate,  the  orientation  of 
the  film  and  all  measurements  are  given  in  units  of  the  STO 
lattice  constants.  In  Figure  2a,  we  show  a  reciprocal  space  map¬ 
ping  of  the  HHL  zone  of  the  film.  The  measurements  were 


performed  at  room  temperature  in  the  reciprocal  spaces  around 
(0.5  0.5  0.5)  where  the  observed  neutron  signal  should  arise 
from  the  magnetic  order  present  in  the  BFO  film,  as  the  STO 
has  no  magnetic  order  and  shows  only  the  nuclear  reflections. 
If  the  spiral  structure  existed  as  observed  in  the  bulk  form,  we 
would  expect  multiple  peaks  in  the  reciprocal  space  around  (0.5 
0.5  0.5)141  as  the  G-type  antiferromagnetism  gives  way  to  a  spiral 
structure  with  propagation  vector  [909]  as  discussed  earlier. 

As  one  can  see  in  Figure  2a,  our  fabricated  (110)  film  clearly 
exhibits  two  reflections  in  the  scattering  zone.  The  spacing 
between  two  peaks  is  0.0053(1)  reciprocal  lattice  units  (rlu) 
which  is  consistent  with  the  bulk  spin  modulation  period  of 
620  A.  Given  that  the  relative  intensities  of  the  two  reflections 
is  similar,  it  is  unlikely  that  they  come  from  two  separate  crys¬ 
tallographic  domains  since  the  X-ray  RSMs  in  Figure  la  and  b 
show  that  this  film,  while  having  two  twins  present,  is  domi¬ 
nated  by  one  twin  domain.  Thus,  we  can  see  that  a  modulated 
structure  is  achieved  in  (1  1  0)  oriented  films.  Importantly, 
these  reflections  are  consistent  with  a  propagation  vector  along 
the  [1  1  -2]  direction  rather  than  the  [1  0  -1]  direction  found  in 
the  bulk.  The  direction  of  the  propagation  vector  was  identified 
from  fits  of  the  observed  intensity  to  two  dimensional  gaussians 
to  determine  the  centers  of  the  reflections.  These  centers  were 
then  used  to  determine  the  propagation  vector  of  [9929]  .  This 
direction  was  also  confirmed  by  measuring 
the  scattering  angle  of  the  observed  reflec¬ 
tions  to  confirm  that  they  were  satellites  of 
the  (0.5  0.5  0.5)  reflection  rather  than  a  (0.5 
0.5  -0.5)  reflection.  Given  the  positions,  they 
are  not  coming  from  possible  other  magnetic 
domains  with  a  propagation  vector  along  the 
<1  1  0>  family  of  directions  being  projected 
into  the  scattering  plane.  Furthermore,  the 
relatively  weak  film  (110)  reflection  was 
found  and  a  measurement  was  performed  in 
that  reference  system  to  verify  the  propaga¬ 
tion  vector.  For  further  testing,  offsets  were 
intentionally  introduced  to  the  instrument 
angles  to  determine  if  instrumental  error 
could  explain  the  position  of  the  reflections. 
All  observations  were  consistent  and  can 
only  be  explained  by  a  propagation  vector  of 
[9929], 

To  examine  possible  formation  of  mag¬ 
netic  domains  in  the  film,  we  show  a  car¬ 
toon  of  the  relative  positions  of  reflections 
in  the  HHL  zone  from  all  possible  magnetic 
domains  with  a  propagation  vector  along  a 
family  of  the  [11-2]  direction  in  Figure  2b. 
These  domains  would  arise  from  rotations 
about  the  three  fold  axis  and  would  result 
in  [1  -2  1]  and  [-2  1  1]  propagation  vectors. 
In  the  HHL  zone  (Figure  2a),  if  magnetic 
domains  are  present,  they  will  be  out  of  the 
scattering  plane;  however,  due  to  the  broad 
vertical  resolution  of  the  triple-axis-spec- 
trometer,  they  would  be  detected  in  the  scat¬ 
tering  plane  as  found  in  previous  neutron 
studies  on  bulk  single  crystals. 1171  However, 
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Figure  2.  Unpolarized  neutron  diffraction  for  the  1  gm-thick  (1 10)  oriented  BiFe03  film.  a,c,d) 
Reciprocal  space  maps  taken  in  HHL,  HKK,  and  HKH  zones,  respectively,  b)  A  cartoon  showing 
the  relative  positions  of  reflections  from  different  magnetic  domains.  Blue  circles  lie  in  the  scat¬ 
tering  plane.  White  circles  represent  the  projection  of  reflections  into  the  HHL  scattering  plane. 
The  cartoon  is  centered  about  a  reflection  of  (0.5  0.5  0.5)  shown  as  a  red  square. 
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in  Figure  2a  there  is  no  indication  of  contribution  from  these 
other  domains  (which  would  be  reflected  in  intensity  present 
at  the  (0.5  0.5  0.5)  ±  (§  §  0)  positions)  implying  that  the  film 
has  a  single  magnetic  domain.  To  more  fully  investigate  this 
point,  we  examined  several  different  scattering  zones.  Repre¬ 
sentative  meshes  are  shown  in  Figure  2c  and  d  for  the  HKH, 
and  HKK  zones,  respectively.  For  both  measured  zones,  a  single 
peak  is  clearly  observed.  If  multiple  domains  were  present, 
instead  of  seeing  a  single  reflection,  we  would  expect  to  see 
multiple  reflections  in  another  zone  in  which  both  satellites 
were  separated.  Thus,  we  conclude  that  this  film  consists  of  a 
single  magnetic  domain. 

Another  reasonable  question  to  ask  is  whether  the  observed 
two  reflections  arise  from  a  modulation  in  the  angle  between 
the  adjacent  magnetic  moments  or  the  amplitude  of  magnetic 
moment,  namely  a  spiral  or  an  amplitude  modulated  magnetic 
structure.  To  address  this,  we  turn  to  polarized  neutron  meas¬ 
urements.  In  polarized  neutron  diffraction,  the  neutron  spin 
state  is  selected  prior  to  scattering  by  the  sample.  After  scat¬ 
tering  from  the  sample,  the  spin  of  the  outgoing  neutron  is 
measured.  We  measure  in  the  channels  in  which  the  neutron 
spin  is  flipped,  or  not  flipped  by  interaction  with  the  sample. 
Neglecting  incoherent  scattering,  the  canonical  equations  gov¬ 
erning  polarized  neutron  diffraction  are  given  byd17^ 
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In  these  equations,  fcf  and  fe;  are  the  outgoing  and  incoming 
wave-vectors  of  the  neutron  respectively.  P  is  the  polarization  of 
the  neutron  and  M  is  the  moment  direction  in  the  sample.  N 
denotes  nuclear  scattering.  Plus  and  minus  indicate  the  polari¬ 
zation  state  of  the  neutron.  Essentially,  these  equations  show 
that  spin  flip  (SF)  scattering,  in  which  the  spin  of  the  neutron  is 
flipped  by  interaction  with  the  sample,  is  sensitive  to  the  com¬ 
ponents  of  the  sample’s  magnetic  moment  perpendicular  to  the 
polarization  of  the  neutron  and  that  non  spin  flip  (NSF)  scat¬ 
tering  is  sensitive  to  the  components  of  the  sample’s  moment 
parallel  to  the  polarization  of  the  neutron.  This  is  schematically 
shown  in  Figure  3. 

For  the  polarized  beam  measurement,  we  focused  on  the  two 
magnetic  reflections  observed  in  the  reciprocal  space  around  (0.5 
0.5  0.5)  in  the  HHL  zone.  Figure  4a  shows  the  scattering  geom¬ 
etry  for  the  measurements.  In  Figure  4b,  we  show  measure¬ 
ments  of  SF  and  NSF  scattering  with  the  neutron  polarization 
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Figure  3.  A  cartoon  of  polarized  neutron  diffraction.  The  neutron  polari¬ 
zation  is  either  in  the  scattering  plane,  or  is  normal  to  it.  The  neutron 
is  sensitive  to  the  component  of  the  magnetic  moment  of  the  sample 
normal  to  the  scattering  wave  vector,  Q  which  lies  in  the  scattering 
plane.  We  denote  this  component  by  M±.  Red  indicates  the  directions 
of  the  electron  moment  to  which  the  neutron  will  be  sensitive,  a)  The 
neutron  polarization,  P,  is  parallel  to  Q  and  we  consider  non  spin  flip 
(NSF)  scattering.  Here,  we  are  not  sensitive  to  the  electron  moment  at 
all  as  discussed  in  the  text,  b)  The  neutron  polarization  is  parallel  to  Q 
and  we  consider  spin  flip  (SF)  scattering.  Here,  the  neutron  is  sensitive 
to  moments  in  the  red  plane,  c)  We  consider  the  case  of  the  neutron 
polarization  perpendicular  to  the  scattering  plane.  In  the  NSF  channel, 
the  cross  section  is  sensitive  to  the  projection  of  the  electron  moment 
along  the  direction  in  which  the  neutron  is  polarized,  shown  in  red. 
d)  In  the  SF  channel,  with  the  neutron  polarization  perpendicular  to  the 
scattering  plane,  the  cross  section  is  sensitive  to  the  projection  of  the 
moment  along  the  red  line  shown  in  the  figure.  We  remind  the  reader  that 
the  grey  plane  shown  is  not  necessarily  the  plane  of  the  film,  but  is  an 
abstract  plane  which  is  normal  to  the  scattering  wave  vector,  Q. 


parallel  to  the  scattering  vector,  Q.  In  this  figure,  we  see  that  SF 
scattering  dominates  NSF  (the  offset  in  backgrounds  is  due  to 
incoherent  scattering  and  the  presence  of  weak  NSF  scattering 
is  due  to  incomplete  polarization  of  the  incident  beam).  This 
tells  us  that  the  magnetic  moment  is  normal  to  the  scattering 
vector  of  (0.5  0.5  0.5)  and  is  thus  either  out  of  the  scattering 
plane,  or  along  the  [1  1  -2]  direction  in  the  scattering  plane, 
although  we  are  not  sensitive  to  components  of  the  moment 
along  the  [111]  direction.  Thus,  Figure  4b  tells  us  that  we  have 
magnetic  scattering,  but  does  not  limit  the  direction.  In  Figure  4c 
we  show  measurements  of  SF  and  NSF  scattering  with  the 
neutron  polarization  out  of  the  scattering  plane.  Again,  we  see 
that  SF  scattering  is  dominant  over  NSF  scattering.  This  leads 
to  the  conclusion  that  the  moment  has  a  component  along  the 
[1  1  -2]  direction,  which  is  one  of  the  axes  in  the  hexagonal  unit 
(with  some  possible  component  along  the  [111]  direction)  cell 
of  BFO.  To  check  this  further,  we  rotated  the  neutron  polariza¬ 
tion  within  the  scattering  plane  and  observed  how  the  intensity 
changes.  In  Figure  4d,  we  show  how  the  intensity  (of  the  right 
peak  from  Figure  4b)  varies  with  the  angle  of  the  sample  guide 
field  rotation.  An  angle  -  0  degrees  corresponds  to  polarization 
along  the  scattering  vector.  The  intensities  have  been  corrected 
for  spatial  inhomogeneity  in  the  spin  transport  using  the  flip¬ 
ping  ratio  of  a  nuclear  peak  measured  over  the  same  range. 
From  this,  we  can  see  that  there  is  variation  and  that  the  NSF 
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Figure  4.  Polarized  diffraction  measurements  of  the  (110)  oriented  BiFe03  film,  a)  Cartoon  of  the  scattering  geometry  for  the  measurements. 
b,c)  Intensity  profiles  taken  in  HHL  zone  with  (b)  P(neutron  polarization)//Q(scattering  vector)  and  (c)  P1.Q.  d)  Changes  in  scattering  intensites  with 
the  polarization  rotated  by  an  angle  0  (measured  in  degrees)  within  the  scattering  plane  relative  to  the  scattering  wave  vector  (that  is,  an  angle  of 
0  corresponds  to  the  neutron  polarization  along  the  scattering  wavevector)  In  (b-d),  the  red  and  black  squares  represent  spin-flip  (SF)  and  non  spin- 
flip  (NSF)  scattering,  respectively.  The  data  have  been  corrected  for  spin  transport  as  described  in  the  text.  Error  bars  where  indicated  are  statistical  in 
origin  and  represent  one  standard  deviation. 


channel  also  shows  variation.  This  is  all  consistent  with  a  mag¬ 
netic  moment  with  a  component  along  the  hexagonal  axes  (the 
[1  1  -2]  direction  in  the  pseudo-cubic  notation)  with  some  com¬ 
ponent  possibly  along  the  [111]  direction.  Thus  we  conclude  that 
in  the  (110)  oriented  him,  we  have  a  single  domain,  modulated 
magnetic  structure  although  we  are  not  able  to  resolve  whether 
the  magnetic  structure  is  a  spiral,  or  a  modulated  structure,  due 
to  our  lack  of  sensitivity  to  a  component  of  the  moment  along 
the  [111]  direction.  However,  we  are  able  to  restrict  the  moment 
to  the  HHL  plane. 

Next,  we  turn  to  the  him  grown  on  the  (111)  oriented  STO 
substrate,  where  a  single  crystallographic  domain  with  rhom- 
bohedral  symmetry  is  maintained!8'  As  with  the  (110)  oriented 
him,  we  start  our  investigations  with  unpolarized  beam  meas¬ 
urements.  Figure  5  shows  a  RSM  around  (0.5  0.5  0.5)  in  HHL 
scattering  zones.  We  see  a  rather  broad  peak,  which  is  broader 
than  the  instrumental  resolution.  Note,  in  contrast  to  the  (110) 
him,  there  are  no  well-separated  reflections  observed  in  con¬ 
trast  to  the  (110)  him  case  in  Figure  2a.  We  also  examined 
several  different  scattering  zones  of  the  him.  Figures  5b  and 
c  show  the  representative  results  for  RSMs  taken  in  HKK  and 
HKH  zones,  respectively.  Clearly  there  is  a  single  broad  peak 
present  in  both  measured  zones.  Thus,  from  the  unpolarized 
beam  measurements,  it  is  hard  to  make  any  conclusions  about 
the  nature  of  the  magnetic  structure.  Thus  again,  we  turn  to 
polarized  beam  measurements. 


This  time,  we  orient  the  him  in  the  (H  K  (H+K)/2)  scattering 
plane,  which  is  dehned  by  the  (1  1  1)  and  (1  -1  0)  reflections  as 
shown  in  Figure  6a.  In  Figure  6b  we  show  SF  and  NSF  meas¬ 
urements  on  the  magnetic  reflection  around  (0.5  0.5  0.5)  with 
the  neutron  polarization  P  parallel  to  the  scattering  vector  Q.  A 
single  peak  is  clearly  observed  for  both  SF  and  NSF  configura¬ 
tions,  indicating  that  there  is  a  (0.5  0.5  0.5)  ordering  wave  vector 
present.  In  Figure  6c,  we  show  SF  and  NSF  measurements  per¬ 
formed  with  P  out  of  the  scattering  plane.  Here,  we  notice  that 
now  NSF  scattering  dominates  SF  scattering.  Furthermore,  we 
note  that  there  are  clearly  two  peaks  present  in  the  SF  channel. 
As  we  only  can  observe  a  separation  of  the  peaks  in  the  SF 
channel  with  P  normal  to  the  scattering  plane,  this  means  that 
the  moment  associated  [0.5  0.5  0.5]  G-type  ordering  wave  vector 
has  a  component  along  the  [1  1  -2]  direction  so  as  to  obscure  the 
separation  seen  in  the  SF  vertical  neutron  polarization  channel. 
As  we  are  insensitive  to  the  component  of  the  moment  along 
the  [111]  direction,  this  places  the  moment  for  this  G-type  order 
in  the  HHL  plane.  This  is  the  reason  for  the  (0.5  0.5  0.5)  broad 
peak  observed  in  the  unpolarized  beam  experiments — there  is 
a  modulated  structure  superimposed  upon  the  magnetic  G-type 
ordering  for  the  (111)  film  (Figure  5).  In  Figure  6d  we  again 
investigated  the  variation  of  the  intensity  with  the  angle  of  the 
in-plane  neutron  polarization.  The  flat  nature  of  the  curve  reveals 
that  the  G-type  ordering  peak  dominates  the  signal  with  its 
moment  having  a  strong  component  along  the  [1  1  -2]  direction. 
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Figure  5.  Unpolarized  neutron  diffraction  for  the  1  gm-thick  (111)  oriented  BiFe03  film.  The  reciprocal  space  mappings  were  taken  in  the  a)  HHL, 
b)  HKK  and  c)  HKH  zones. 


Unfortunately,  we  can  say  relatively  little  about  the  modu¬ 
lated  structure  for  the  (111)  film.  However,  one  thing  which  is 
clear  from  the  observations  in  Figure  6c  is  that  the  modulated 
magnetic  moment  is  comparable  in  magnitude  to  the  G-type 
ordering,  and  has  a  strong  component  along  [1  -1  0],  lying 
within  the  plane  containing  [1  -1  0]  and  [1  1  1],  Finally,  we  note 
that  despite  looking  for  evidence  of  chirality  by  distinguishing 
the  -I—  channel  from  the  — f  channel,  we  could  not  observe  any 
asymmetry.  However,  this  may  simply  be  due  to  the  presence  of 
chirality  domains  in  the  film  that  were  not  present  in  the  bulk. 


Our  observations  have  demonstrated  that  it  is  possible  to 
recover  a  modulated  structure  in  the  thin  films.  Examples  of 
possible  modulated  magnetic  structures  consistent  with  our 
observations  are  shown  in  Figure  7.  In  addition,  we  have  found 
that  the  substrate  orientation,  which  results  in  different  strain 
relaxation  processes  leading  to  the  different  thin-film  structural 
phased  has  a  bearing  on  the  magnetic  structure  in  the  film. 
Our  preliminary  neutron  diffraction  studies  for  the  (001)  ori¬ 
ented  film  shows  that  there  is  only  a  single  reflection  around 
(0.5  0.5  0.5),  revealing  that  the  magnetic  structure  of  the  (001) 


Figure  6.  Polarized  diffraction  measurements  of  the  (111)  oriented  BiFe03  film,  a)  Cartoon  of  the  scattering  geometry  for  the  measurements. 
b,c)  Intensity  profiles  taken  with  (b)  P(neutron  polarization)//Q(scattering  vector)  and  (c)  P±Q.  d)  Changes  in  scattering  intensities  with  the  polariza¬ 
tion  rotated  by  an  angle  9  (measured  in  degrees)  within  the  scattering  plane  relative  to  the  scattering  wave  vector  (that  is,  an  angle  of  0  corresponds 
to  the  neutron  polarization  along  the  scattering  wavevector)  In  (b-d),  the  red  and  black  squares  represent  spin-flip  (SF)  and  non  spin-flip  (NSF)  scat¬ 
tering,  respectively.  The  data  have  been  corrected  for  spin  transport  as  described  in  the  text. 
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the  same  mechanism  for  device  applications.  Even  if  the  mag¬ 
netic  structure  is  amplitude-modulated,  it  still  opens  the  pos¬ 
sibility  of  novel  device  applications.  The  electric  field  induced 
magnetization  switching  of  exchanged  coupled  Co  layer  on 
a  (001)  oriented  BFO  film  previously  reported'19!  points  to  a 
mechanism  other  than  one  involving  the  spiral.  The  implica¬ 
tion  of  this  to  the  present  observation  is  that  the  ferroelectric 
domains  in  BFO  films  could  couple  to  the  amplitude-modulated 
antiferromagnetic  domains  and  that  this  coupling  mechanism 
can  be  utilized  for  thin  film  device  applications. 


(b)  (1 1 1)  oriented  film 


Figure  7.  Cartoon  of  possible  magnetic  structures  for  (110)  and  (111) 
oriented  BiFe03  thin  films,  a)  A  possible  magnetic  structure  for  (110)  film 
where  the  magnetic  moments  lies  in  the  HHL  plane  containing  the  [111] 
and  [11-2]  directions.  The  propagation  vector  is  along  the  [11-2]  direc¬ 
tion.  b)  A  possible  magnetic  structure  for  the  (111)  film  where  the  G-type 
ordering  occurs  in  the  plane  containing  the  [111]  and  [1 1-2]  directions, 
and  the  modulated  components  appears  along  the  [-110]  direction  with 
the  propagation  vector  is  along  the  [11-2]  direction. 


3.  Conclusions 

We  have  shown  that  is  possible  to  recover  a  modulated  mag¬ 
netic  structure  in  epitaxial  BiFeOj  thin  films.  Furthermore, 
we  have  found  that  the  magnetic  structure  strongly  depends 
on  the  substrate  orientation  of  the  film.  For  the  (001)  oriented 
substrate,  we  found  that  there  was  only  commensurate  G-type 
antiferromagnetism  present  even  for  films  as  thick  as  1  pm. 
Whereas,  1  pm-thick  films  grown  on  (110)  and  (111)  oriented 
STO  substrates  evinced  modulated  magnetic  structures.  In  the 
(110)  film,  we  found  a  single  magnetic  domain  with  a  modu¬ 
lated  structure,  whereas  both  G-type  and  modulated  magnetic 
structures  were  found  in  the  (111)  oriented  film.  The  recovery 
of  a  modulated  magnetic  structure  in  thin  films  opens  the 
door  to  novel  coupling  between  BFO  films  and  magnetic  layers 
deposited  on  top  of  them. 


oriented  film  is  G-type  antiferromagnetic  structure  (see  Sup¬ 
porting  Information).  However,  for  the  (110)  oriented  film,  a 
single  domain,  modulated  magnetic  structure  is  observed  in 
which  spins  lie  in  the  HHL  plane.  The  propagation  vector  asso¬ 
ciated  with  this  order  is  different  from  that  found  in  the  bulk. 
For  the  (111)  oriented  film,  a  modulated  structure  is  observed, 
but  the  G-type  order  is  also  present.  It  is  interesting  to  point  out 
that  the  structural  symmetry  for  the  (001)  and  (110)  oriented 
film  has  been  lowered  to  the  monoclinic  Ma  and  M],  phase, 
respectively^8'.  This  symmetry-lowered  structural  phase  could 
lift  the  degeneracy  in  the  easy  magnetization  axis  and  change 
the  ground  state  for  the  magnetic  structure.  This  could  be  a 
reason  for  the  single-domain  modulated  magnetic  structure 
observed  in  the  (110)  film.  Recent  first-principles  calculations'3! 
also  suggest  that  the  monoclinic  distortion  may  be  the  reason 
for  creating  the  easy  magnetization  axis  in  BFO.  This  also  sup¬ 
ports  our  experimental  observations. 

Another  implication  of  our  observation  is  the  possibility  of 
growing  magnetoelectric  devices  in  film  form  based  on  the 
spiral  magnetic  structure.  The  recent  work  of  Lebeugle  et  al.'18' 
has  shown  that  it  is  possible  to  electrically  switch  the  easy  mag¬ 
netization  axes  of  permalloy  deposited  on  single  crystals  of 
BFO.  They  suggest  that  the  mechanism  is  strongly  related  to 
the  presence  of  the  spiral.  We  have  shown  that  changing  the 
film  growth  orientation  is  sufficient  to  re-establish  a  modulated 
magnetic  structure  in  BFO  thin  films.  Further  work  needs  to  be 
done  to  determine  whether  the  modulation  comes  from  a  spiral 
or  an  amplitude  modulated  structure.  If  it  turns  out  that  the 
structure  is  a  spiral  as  in  the  bulk,  we  should  be  able  to  exploit 


4.  Experimental  Section 

BiFe03  thin  films  with  thicknesses  up  to  1  pm  were  grown  epitaxially 
on  (001),  (110),  and  (111)  oriented  SrRu03-buffered  SrTi03  substrates 
by  pulsed  laser  deposition.  During  the  deposition,  the  substrate 
temperature  and  oxygen  pressure  were  kept  at  590  °C  and  25  mTorr 
(3.33  Pa).  For  structural  characterization,  a  conventional  4-circle  XRD 
(Bruker  D8)l20!  was  used.  Ferroelectric  properties  were  characterized  with 
a  Radiant  Premiere  III20!  loop  analyzer  at  room  temperature. 

Unpolarized  and  polarized  neutron  diffraction  measurements 
were  performed  respectively  on  the  BT9  and  BT7  thermal  triple-axis 
instruments  at  the  NIST  Center  for  Neutron  Research.  The  measurements 
performed  on  BT9  used  relatively  tight  Soller  collimations  of  40'-lT- 
S-40'-120'  full-width-at-half  maximum  (FWHM).  We  placed  PG  filters 
both  before  and  after  the  sample  to  suppress  higher  order  wavelengths. 
For  the  polarized  beam  measurements  performed  on  BT7,  we  used 
slightly  more  relaxed  Soller  collimations  of  50'-25'-S-50'-open  FWF1M. 
Polarization  was  achieved  through  the  use  of  3Ple  cellsl21!  placed  before 
and  after  the  sample.  Flipping  ratios  as  high  as  10  were  achieved  initially. 
However,  as  the  polarizion  of  the  cell  decays  with  time,  flipping  ratios 
lowered  with  time.  Cells  were  replaced  during  the  experiment  to  maintain 
acceptible  flipping  ratios.  Measurements  presented  in  the  figures 
are  close  enough  in  time  that  corrections  for  the  lowered  polarization 
of  the  cell  are  unnecesary.  All  measurements  were  done  at  room 
temperature.  Measurements  on  a  BFO  film  grown  on  (0  0  l)-oriented 
STO  were  also  performed  on  the  DIO  four-circle  diffractometer  at  the 
ILL  (Grenoble,  France),  in  a  beam  of  wavelength  2.359  A  from  a  PG 
[0  0  2]  monochromator  on  a  thermal  neutron  guide.  A  PG  filter  was 
placed  before  the  sample  to  suppress  higher-order  contaminations. 
The  resolution  at  focusing  was  20’.  Data  were  collected  with  a  position 
sensitive  detector  with  pixel  spacing  equivalent  to  12’  in  both  directions, 
but  the  true  resolution  was  dominated  by  the  dimensions  ofthe  substrate 
face  and  the  1 .7°  vertical  focusing  of  the  monochromator. 
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